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MICROSTABILITY  OF  A FOCUSSED  ION  BEAM  PROPAGATING 
THROUGH  A Z-PINCH  PLASMA 


I.  INTRODUCTION 

In  order  to  utilize  intense  ion  beams  in  a pellet  fusion  device, 
it  is  generally  considered  necessary  to  propagate  the  beam  a distance 
L (about  10  m)  from  the  acceleration  region  while  employing  some 
focussing  scheme  to  deliver  the  beam  to  the  target.  For  heavy  ion 
beams  (e.g.  ^0  GeV  J+  ions),  field -free  ballistic  focussing1’ 2 to  the 
target  has  been  proposed.  For  light  ion  beams  (e.g.  10  MeV  protons), 
one  concept3  involves  focussing  the  beam  just  outside  the  diode  region 
by  known  techniques4  and  then  injecting  the  focussed  beam  into  a z- 
pinch  plasma  channel.  The  beam  remains  focussed  during  propagation  by 
the  magnetic  field  in  the  plasma  channel  and  is  delivered  on  target 
with  no  additional  focussing.  This  scheme  for  light  ion  beams  is 
considered  in  this  paper.  Since  the  plasma  channel  must  remain  intact 
during  transit  of  the  beam,  a knowledge  of  the  dependence  of  the  linear 
growth  rates  of  unstable  perturbations  on  the  parameters  of  the  system 
is  essential.  Here  the  concern  will  be  with  microinstabilities. 

It  will  be  assumed  that  the  z -pinch  channel  is  produced  in  a MHD 

stable  configuration.  The  question  of  the  effects  that  the  passage  of 

the  beam  may  have  on  the  MHD  equilibrium  and  stability  of  the  beam- 
Note:  Manuscript  submitted  May  26,  1978. 


plasma  system  is  an  important  one.  For  example,  the  primary  and  return 
current  heating  of  the  plasma  by  the  beam  may  significantly  modify  the 
pressure  profile  of  the  background  plasma  and  thus  upset  the  MHD  equi- 
librium. This  question  is  not  addressed  here,  but  instead  it  will  be 
assumed  that  the  plasma  channel  remains  in  a stable  MHD  equilibrium 
during  the  passage  of  the  beam  in  order  to  study  the  microinstabilities 
of  the  system. 

Because  the  ion  beam  streams  through  the  background  plasma,  the 
system  may  be  subject  to  velocity -space  instabilities.  Two  modes  will 
be  investigated.  One  mode,  which  will  be  referred  to  as  the  e-b  mode, 
involves  the  interaction  of  beam  ions  and  background  electrons  and  the 
other  mode,  the  e-i  mode,  involves  the  interaction  of  background  elec- 
trons and  background  ions.  Relative  motion  between  the  stationary  back- 
ground ions  and  the  background  electrons  arises  from  both  the  z-pinch 
electron  current  and  the  electron  return  current  driven  by  the  ion 
beam;  however,  the  return  current  provides  the  larger  contribution 
since  the  ion-beam  current  greatly  exceeds  that  establishing  the 
channel. 

An  electrostatic  stability  analysis  is  used  for  simplicity  and  the 
wavevector  k is  chosen  to  be  parallel  to  the  direction  of  beam  propaga- 
tion (i.e.  k = ke  ) . If  either  of  the  two  modes  that  are  considered  is 
— z 

sufficiently  unstable  to  create  strong  microturbulence,  the  resulting 
anamolous  resistivity  could  seriously  affect  the  current  neutralization 
of  the  beam  which,  in  turn,  might  disrupt  beam  transport  because  of  the 


perturbed  magnetic  guide  field.  Thus  it  is  important  to  know  in  which 
parameter  regimes  these  modes  are  stable  in  order  to  design  a system 
which  will  maintain  good  current  neutralization.  It  should  be  noted 
that  the  high  density  of  the  z-pinch  plasma  required  to  insure  good 
beam  neutralization5  (n  =“  5 x 10la  cm  3)  implies  that  collisional 
processes  will  play  an  important  role  in  the  stability  analyses. 

Two  models  will  be  considered  for  the  z-pinch  plasma  channel.  In 
Section  II,  a surface  current  model  will  be  used  to  describe  the  plasma 
channel  and  in  Section  III  a uniform  current  model  will  be  used. 
Although  a uniform  current  model  more  closely  approximates  the  z-pinch 
plasma  which  is  envisioned  in  this  focussed  ion  beam  propagation  scheme, 
the  surface  current  model  is  treated  here  as  well  for  the  sake  of 
clarity.  In  both  cases,  expressions  for  the  linear  growth  rates  of 
the  e-b  mode  and  the  e-i  mode  will  be  derived.  From  these  results,  it 
will  be  demonstrated  that  it  is  possible  to  propagate  a focussed  ion 
beam  without  the  development  of  microinstabilities.  This  conclusion 
and  other  results  will  be  discussed  in  Section  IV. 

II.  MICROSTABILITY  ANALYSIS  - SURFACE  CURRENT  MODEL 

The  case  where  the  z-pinch  current  flows  on  the  surface  of  the 
plasma  column  will  be  considered  first.  For  simplicity,  no  applied 
axial  magnetic  field  will  be  considered.  Thus  the  beam  propagates  in 
a field-free  region  and  beam  ions  only  encounter  the  azimuthal  field, 

Ba,  when  they  are  at  the  edge  of  the  beam.  Writing  the  field  as 

0 r < a! 

Bpa/r  r > aj  , (1) 


3 


where  B is  the  magnitude  of  the  peak  magnetic  field,  the  vector 
potential,  A , becomes 


A = e 


-B^a^n  (r/a) 


r < a 


r > a 


The  beam  ions  have  straight  line  orbits  in  the  interior  of  the  beam, 
and  are  reflected  back  into  the  interior  when  they  encounter  the 
azimuthal  field  at  the  edge  of  the  beam.  If  the  beam  distribution 


function  is  written  as 


where  Vz  is  a constant  streaming  velocity  of  the  beam, 

= n^(vz  + eAz/m^c)  is  the  social  canonical  momentum  and  is  the 
beam  thermal  energy,  then  from  Eq,  (2) 


r < a 


(|)S  r > a 


where  is  the  beam  density  in  the  interior  of  the  beam, 

0 * 2vzu,Cba/u§>  5 eB^/n^c  and  = m^u^/2.  The  condition  for  the 
beam  to  be  well  confined  by  B.  (i.e.  the  beam  density  falls  off  rapidly 
for  r > a)  is  just 

0 > 1 . 


(5) 


The  same  condition  on  0 is  obtained  from  the  beam-ion  orbit  equations 


i 


Li 


by  requiring  that  (rffl  - a)/a  < 1,  where  rm  is  the  maximum  radial  excur- 
sion for  any  given  beam  ion. 

For  a 10  MeV  proton  beam  (V  = 4.1+  x 109  cm/sec)  with  a thermal 

z 

energy  of  T^  = 0.1  MeV  (u^  = 4.1+  x 10s  cm/sec),  Eq.  (5)  shows  that 

B > 4.2  kG  is  required  in  order  to  propagate  a focussed  beam  of 
P 

radius  a = 0.5  cm.  Typically,  B =“  40  kG  will  be  used  (i.e.  0 =“  10). 

A typical  system  will  also  have  a beam  density  of  n^  < 1 x 1016  cm"3 

(i.e.  n^/n^  <<:  !)•  Thus,  with  full  current  neutralization  the  electrons 

drift  relative  to  the  stationary  ions  with  an  axial  velocity  of 

V “■  n, V /n  =-  9 X 10®  cm/ sec. 
e o z p 

The  conservation  of  axial  canonical  momentum  shows  that  a beam 
ion  moves  with  a larger  axial  speed  when  r > a than  when  r < a.  In 
fact,  for  r > a 


Avz  = vz  "Vz  = tweba^l(r/a)* 


(6) 


During  the  time  it  takes  for  a beam  ion  to  reflect  off  the  magnetic 

field  barrier  (r  > a),  this  increased  axial  speed  moves  the  beam  ion  a 

distance  AZ  =*■  4au^0/Vz  ahead  of  the  position  it  would  have  if  it  were 

streaming  in  the  interior  of  the  beam  (r  < a)  with  a speed  V . If 

z 

kAz  > 1,  then  this  process  is  equivalent  to  an  effective  collisional 

process  knocking  a beam  ion  out  of  phase  with  a wave  moving  with  a 

phase  velocity  tu/k  *"  V each  time  the  beam  ion  reflects  off  the 

z 

magnetic  field.  This  process  is  not  important  here  because  the  effec- 
tive collision  time,  t *■  a/u^,  is  much  longer  than  the  time  scales  of 


5 


interest  for  the  e-b  mode,  however,  a similar  process  will  be  important 
when  considering  the  e-b  mode  for  the  uniform  current  case  discussed 
in  Section  III,  since  in  that  case  the  beam  ion  motion  is  influenced 
by  Bq  at  all  r,  not  just  r > a. 

Considering  the  above  results,  a beam-plasma  system  with  the 
Z-pinch  current  flawing  on  the  surface  of  the  plasma  channel  can  be 
assumed  to  be  infinite  and  uniform  for  the  purposes  of  this  stability 
analysis.  This  is  also  a reasonable  assumption  since  it  can  be  shown 
that  ka  » 1 in  the  parameter  regimes  of  interest.  The  electrons  are 
treated  as  a cold  collisional  fluid  and  the  beam  ions  are  treated 
kinetically.  The  appropriate  dispersion  equation  for  the  e-b  mode  is 
(see  Appendix  A) 


1 


or 


Pe 

uu(u)+iv  ) 


- h + 
k2uf  \ 


i/TT(u)-kV  ) 
z 

K 


exp 


(T) 


-a  , 3/a 

where  = 1.5  x 10  ne^ei/e  ; electron  collision  frequency, 

\ . is  the  Coulomb  logrithm  for  electron-ion  collisions  and  u>  is  the 
ei  per 

plasma  frequency  for  species  a (T  is  expressed  in  eV).  Solving  for 
the  maximum  growth  rate  uu®,  of  the  e-b  mode,  one  obtains 


uu 


r 


(9) 


(9) 


v i—  uu2,  V2 


pe  b 


where  cu  = iuK.  The  condition  for  stability  of  the  e-b  mode  is  just 


[± MVWVf  >:l 
s «*.  VV 


For  the  parameters  user1  previously,  the  e-b  mode  is  stable  for 
Tg  ^ 50  eV.  Since  it  is  expected6  that  the  primary  and  return  current 
heating  will  raise  the  electron  temperature  no  higher  than  this  value, 
the  e-b  mode  is  stable  for  the  case  of  surface  currents. 

Because  good  current  neutralization  requires  n^/n^  « 1,  the 
streaming  velocity,  V , of  the  electrons  relative  to  the  stationary 
background  ions  is  actually  subthermal  (V  < u , where  T = m u2/2) 

S 6 S S S 

for  T ^ 1 eV.  Thus,  the  electrons  can  be  considered  as  a warm  colli - 
e 

sional  fluid  when  studying  the  e-i  mode.  Since  T > T . fu,  « u ), 

e iv  i e ’ 

the  ions  can  be  treated  as  a warm,  collisional  fluid  as  well.  In  the 

momentum  transport  equations,  the  collisional  drag  force  on  species  a 

is  expressed  as  -m  n v _(v  - v „ ) and  the  thermal  pressure  gradient 
or  or  ~0 

term  is  expressed  as  Thus  the  dispersion  equation  for  the  e-i 

mode  becomes 


^pjC^d  " k2uf/2)+  ^pg ( aj2  ~ 

(i»u(i[vev^  + (oj  + iv^  - k2u2/2uu)(uj(i  + ivg  - k2u2/2u>g)] 

where  uk  s a>  - kV  , v.  = m v ,/m.  and  v was  defined  earlier  (see 
d e 7 i e e 1 e 

Appendix  A for  details  of  the  derivation). 
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Solving  Eq.  (12)  for  the  real  and  imaginary  parts  of  u),  one 


obtains  for  T.  < T 
i e 


(JO 


= kC./(l  + k2X|)^ 


r s 


and 


u>,  = 


1 2(1  + k2X2)  L 


1 + 


Ve/Cs  - 1/(1  + k^2)^ 


(1  + k2!2)* 


(13: 


(IV 


where  C2  = T /m.  is  the  ion  sound  speed  and  Xf  = ua/2uo2  is  the 
s e'  i d e'  pe 

electron  Debye  length.  In  deriving  this  result  it  was  assumed  that 

k^2  > v UK  > jo?  and  oo  > ku.  > v . Since 
e e d d i i 


o»i  < 0 for  all  k (15 

the  e-i  mode  is  stable. 

Sumnarizing  the  results  obtained  in  this  section,  it  was  found 
that  e-i  mode  is  stabilized  by  collisional  damping  and  that  the 
e-b  mode  is  stable  because  the  electrons  are  not  heated  above  a 
critical  temperature  determined  by  Eq.  (ll)  for  the  parameters  of 
interest.  These  results  apply  to  the  case  where  the  z-pinch  current 
flows  on  the  surface  of  the  plasma  column.  If  the  z-pinch  current 
is  distributed  uniformly  across  the  plasma  channel  the  azimuthal 
magnetic  field  within  the  channel  can  have  an  additional  stabilizing 


influence . 


in.  MICROSTABILITY  ANALYSIS  - UNIFORM  CURRENT  CASE 


If  the  z-pinch  current  is  uniformly  distributed,  then  the  azimuth- 
al magnetic  field  increases  linearly  with  radius  inside  the  beam.  The 
beam  ion  orbits  are  now  nearly  sinusoidal  when  the  ion  enters  the 
channel  at  a small  angle  to  the  plasma  axis.  For  mathematical  con- 
venience, a slab  model  will  be  used  here.  This  is  appropriate  for 
beam  ions  with  small  angular  momentum  since  motion  then  occurs  in  a 
plane.  If  the  beam  is  streaming  in  the  z direction,  then  the  aximuth- 
al  field  will  point  in  the  y direction  and  can  be  written  as 

B^/a  - a < x < a 

Bp  (x|  > a 

The  vector  potential.  A,  in  this  geometry  becomes 

I-B  x2/ 2a  - a < x < a 

P 

-Bp (x -a/2)  j x | > a 

(17) 


Again  writing  the  distribution  function  as  in  Eq.  (3),  f^  can  be  put 
in  the  following  form 


J 


where  again  0 = ^Vz<JUcba/ub’  aucti  = eB^/m^C  and  again  the  condition  for 

confinement  of  the  beam  is  just  0/2  ^ 1.  The  beam  ion  orbits  are 

affected  by  the  axial  electric  field  (|e  | = m v v/e)  which  is  re- 

quired  to  drive  the  z-pinch  current,  however,  on  time  scales  shorter 

than  the  beam  pulse  length  (t  < 50  nsec)  the  effect  of  E is  negligi- 

z 

ble.  Using  Eq.  (1 6)  for  B , the  beam  ion  equation  of  motion  can  be 

•7 

solved  approximately  in  terms  of  trigonometric  functions. 


x'  = x cosu)  t + — sinou  t 
0 U>  o 

o 


(19) 


k . ; 


y'  - y + v t 


U) 

+ BaT 

z 


/ V2 

lx2  - -| 

\ ^o 


(20) 


(JO2  /V2  V* 

z'  =z+LVz  + ctb'xT 


vf  \ XV 

Jsin  2'juqt  - (cos2tuor-l),  (21) 


where  uj  = (uj  , V /a)^  and  (x1 , v' ) = (x,  v)  at  t = t'  - t = 0. 

^ CD  Z 


This 


solution  is  obtained  by  replacing  v'  with  v in  the  x -component  of  the 

z z 


equations  of  motion  and  is  therefore  applicable  when  '^yJv'z  « 1»  In 
this  case,  an  expressior 
obtained  frcm  Eq.  (2*0, 


this  case,  an  expression  for  Av  (averaged  over  a 2iu  period)  can  be 

z o 


Av,  = f 


V'  - V 


uj  , a 
cb 


(25) 


where,  for  example,  an  ion  starting  at  t = 0 from  the  point  of  maximum 
excursion  (x'  = x = a and  v ' = v =0)  will  have  an  axial  velocity 

X X 

given  by  v'(t)  = v + (a2u)2Av  )(cos2u)  t-1).  Using  the  parameters 
z z o z o 

from  the  example  in  Section  II,  as  well  as  Eq.  (9)  for  km  derived  for 
the  e-b  mode,  it  is  found  that  kmAv  » uu?  , so  that  the  average  beam 

Z 1 

ion  moves  out  of  phase  with  the  wave  even  before  one  e -folding  time. 

In  addition,  uu  L/v  » 1,  so  that  the  beam  ions  complete  many  oscilla- 
o z 

tions  while  propagating  down  the  plasma  channel.  Because  of  the 
oscillatory  motion  in  the  z -direction  the  beam  ions  are  constantly 
moving  in  and  out  of  phase  with  waves  moving  with  phase  velocity 
iur/k  -•  V . This  process  is  responsible  for  a reduction  in  the  growth 
rate  of  the  e-b  mode,  which  will  now  be  derived. 

Using  the  method  of  characteristics7,  the  perturbation  to  the 
beam  distribution  function  can  be  written  as 

*«.(£»£»*)  = J (tf-  -^r  + 1^)  exp[i(kZ-o,T)]dT  , (26) 
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for  electrostatic  perturbations  of  the  form  0 = 0(x)exp[i(kz-iut)] . 

Here  Z = z ' - z , k = 0 and  the  x dependence  has  not  been  Fourier 

y 

analyzed.  Since  the  solution  of  the  eigenfunction  equation  for  0(x) 
is  difficult  and  of  little  present  value,  the  solution  will  be  approxi- 
mated by 


x2 

1 - — 


< X < 


0(x)  = 0 


so  that  the  perturbation  field  is  confined  within  the  beam.  Here 
s (2/0 )^a  is  the  beam  radius  and  a^  < a for  0 > 2.  The  perturba- 
tion is  chosen  to  peak  on  axis  in  order  to  model  the  zero -order 
eigenmode  (no  nodes).  The  expression  for  f^,  Eq.  (26),  is  now 
integrated  over  velocity  space  to  obtain  an  expression  for  the  beam 
density  perturbation,  n^,  which  is  inserted  into  Poissons  equation 


-7a0  = 1+rre(nbl  - ngl) 


along  with  the  appropriate  expression  for  n^.  Again  the  electrons 
are  taken  to  be  a cold  collisional  fluid  when  treating  the  e-b  mode. 

Taking  a weighted  spatial  average  of  Poisson's  equai  Ion  results 
in  an  approximate  algebraic  dispersion  equation 


a)(<u+iv 


Mm-kVz 


▼to 


The  details  of  the  calculations  leading  to  this  dispersion  equation 
are  given  in  Appendix  B. 


Eq.  (29)  has  the  same  form  as  Eq.  (7)  except  that  u)^  is  replaced 
ujp^Rb  where  R^  reflects  the  reduction  of  beam  ions  in  resonance  with 
the  wave  as  discussed  earlier. 

\ - VoVK  • (30> 


1 


I 

I I 


Note  that  is  the  ratio  of  the  wavelength  of  the  perturbation  to  the 
amplitude  of  the  oscillatory  motion  in  the  z direction  of  the  beam 
ions  (see  Eq.  (21)  for  the  oscillations  in  z at  2u)q).  Also  note  that 
for  the  collisional  electrons  |uj(ux-ive)|  » |eBp/mec|2. 

The  condition  for  stability  of  the  e-b  mode  for  the  uniform 
current  case  becomes 


4 

In  this  case  the  e-b  mode  is  stable  for  Tg  <1.0  keV  when  the  parame- 
ters used  in  the  example  in  Section  II  are  inserted  into  Eq.  (31). 

This  is  well  above  the  temperatures  that  are  expected.  Thus  the  uni- 
form current  case  is  more  desirable  than  the  surface  current  case  with 
respect  to  the  stability  of  the  e-b  mode. 

The  analysis  of  the  e-i  mode  for  the  uniform  current  case  differs 
very  little  from  the  analysis  for  the  surface  current  case  since 
uj  « ku  , <ju,  v and  oi  . « <u,  v. . Although  the  collisional  fluid 
treatment  of  the  electrons  and  ions  do  not  include  finite  Larmor 
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radius  effects,  such  effects  will  only  further  stabilize  the  system. 
Thus  the  e-i  mode  will  be  stable  for  the  uniform  current  case  just  as 
in  the  surface  current  case. 

IV.  CONCLUSIONS 

The  microstability  analysis  presented  here  treated  the  e-b  and 
e-i  modes  of  a focussed  ion  beam-plasma  system.  The  high  density 
(n^  ~ 5 x 1018  cm-3)  required  for  the  background  z-pinch  plasma  pro- 
vides for  a high  frequency  of  collisions  between  the  background 
electrons  and  ions.  It  was  found  that  this  results  in  the  stabiliza- 
tion of  the  e-b  mode  by  collisional  damping,  even  for  reasonably  high 
electron  temperatures.  In  the  example  presented  in  Section  H,  the 

e-b  mode  was  stable  for  T £ 50  eV  for  the  surface  current  case. 

e 

When  the  z-pinch  current  is  distributed  uniformly  over  the  plasma 
channel,  it  was  found  that  the  e-b  mode  was  even  further  stabilized  by 
the  oscillatory  motion  (in  the  direction  of  beam  propagation)  exhib- 
ited by  the  beam  ions.  In  this  case  the  e-b  mode  was  stable  for 
Tfi  ^ 1.0  keV  for  the  same  parameter  as  used  in  the  previous  example. 
Since  the  electrons  are  expected  to  be  heated  no  higher  than  50  eV, 
z-pinch  current  flow  in  the  plasma  channel  ensures  stability  of  the 
e-b  mode.  Uniform  current  flow  is  anticipated,  because  the  magnetic 
diffusion  time  in  a current  channel  1 cm  in  diameter  is  much  less  than 
the  many  nsec  required  to  establish  the  channel.  The  low  channel 
currents  required  to  establish  the  kO  kG  azimuthal  fields  are  suf- 
ficient to  heat  the  channel  to  more  than  a few  eV  before  the  ion 
beam  enters  the  channel. 


Electron-ion  collisions  are  also  responsible  for  stabilization  of 

the  e-i  mode.  The  condition  n^/ « 1,  required  for  good  charge  and 

current  neutralization  of  the  beam,  implies  the  Vg  < ug  for  a typical 

system.  This  is  important  since  the  e-i  mode  can  be  unstable  if 

V > u . The  results  for  the  e-i  mode  are  basically  unchanged  if  the 
e e 

current  flows  in  the  interior  of  the  beam  channel. 

Thus  both  the  e-b  and  e-i  mode  are  stable  for  a typical  system  if 
the  focussed  ion  beam  is  propagated  through  a z -pinch  plasma  channel 
where  the  z -pinch  current  flows  in  the  interior  of  the  channel. 


APPENDIX  A 


I 


li 


Treating  the  plasma  electrons  and.  ions  as  collisional  fluids,  the 
linearized  fluid  equations  describing  their  motion  under  the  influence 
of  an  electrostatic  perturbation  of  the  form  exp  i(kz-ujt)  become 


(u>-kV  )n  = n k*v 
or  on  or -cr  i 


(A.l) 


and 


m (oo-kV  )v  = q kn0,  + T kn  /n  - iv  m (v  - vQ),  , (A. 2) 

o'  a'—ai  1 o~  or  o o o ~o  i ’ ' 


where  V ^ is  the  streaming  velocity,  is  the  thermal  energy  and  is 
the  collision  frequency  for  species  o(o  = e,i),  and  the  subscript  ”i" 
signifies  a perturbation  quantity.  Solving  for  n in  terms  of  0* 


n q k20,/®  (uj  - kV  ) 
o o v o'  cr 


oi 


+ j^uj  + iv^  - k^2/ 2ir|  £iu  - kVg  + ivg  - )hsm^/2{w  - kVg  )J^ 


, (A.3) 


where  = 0 and  = m^u^/2.  This  expression  for  n , Eq.  (A.3),  is 
inserted  into  Poisson's  equation  in  order  to  obtain  the  dispersion 
equation. 

For  the  e-b  mode  the  beam  ions  are  treated  kinetically  so  that 


11  -b  J 


k • df^/dv 
uj-k-v 


0 d3v 

l — 


(A.4) 
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where  f^  is  given  in  Eq.  (M.  Rewriting  Eq.  (A.1*)  in  terms  of  the 
plasma  dispersion  function,  Z(|),  results  in 


V 


[1  + 5Z(5)]  , (■ 

Vb 


where  § = (co-kV^/ku^.  For  a warm  beam  § < 1 and  the  power  series 
expansion  of  Z(§)  is  appropriate8. 


APPENDIX  B 


The  ion  beam  contribution  to  Poisson's  equation,  Eq.  (28),  is 


given  by 


d3v  £ 


(B.l) 


where  is  given  in  Eq.  (26).  Taking  the  weighted,  spatial  average  of 
Poisson's  equation,  the  desired  expression  for  the  dispersion  equation, 


Eq.  (29),  becomes 


<v> 


(B.2) 


Assuming  8/2  > 1,  so  that  for  all  x 


♦(\-V 


(B.3) 


Eq.  (B.2)  becomes 


<nbi>= 


:(1’7r)fi‘Vz)  exp{'^;  vx+vy  + (\-Vz)S  » 


(B.4) 


where  Z = z'  - z , r = t ' - 1,  and  x'and  v’  are  given  in  Eqs.  (19) -(24). 
Then  using  the  Bessel  function  indentities9 


g±ia  sin  bt_  ^m(a)cos(2mbt)  ± 21^  J2m+1(a)  sin[ (2m+l)bt] 


and 


gia  cos  bt  _ + 2 i j(a)cos(mbt) 


(B.5) 

(B.6) 


to  simplify  exp[ikZ],  one  obtains  a sun  of  integrals,  where  each 
integral  in  the  sum  has  the. form 


I(p>s  ;T),u,ct) 


where  = 8uu  V ./k  and  where  the  integration  over  x has  been  extended 
to  infinity  since  exp(-0x2/ 2a2)  -*  0 for  |x|  > a^  = (2/6 )^a  when 
d/2  > 1;  also  p,s,71,n  and  a are  all  integers.  Because  of  the  oscilla- 
tory nature  of  the  Bessel  functions  and  the  properties  of  the 
exponential  function,  the  major  contribution  to  I occurs  at  x « 0 and 
vx  « 0.  In  fact,  l(o,o;Tl,n,a)  » l(p,s ;?l>^,ff)  for  any  p > 0 and  s > 0 

Furthermore,  the  concern  is  with  the  resonance  at  uo-kV  ±2auu  where 

z o 

I 

m = 0.  The  largest  contribution  at  this  resonance  is  obtained  when 
11  = n = a = 0.  Performing  the  integration 


I(o,o;o,o,o)  = Jn 


(B.8) 


where 


-CD  < 

/d3/ 

O " c\ 


dt  (exp 


['1 

(£*•■ ) 

J0(s2)JQ(t2)  JQ(2st )cos(2st) 

L % 

' 'J 

(B.9) 


For  typical  values  of  8 and  V (9  = 10.0  and  V^/u^  «-  1.2  x 10-2 
for  a^  = 0.5  cm,  °*  ko  kG  and  all  other  parameters  the  same  as  used 
previously),  crQ  = O.kk. 

After  performing  the  remaining  t,  v and  v integrations  in  that 

y z 

order,  one  finds 


Une  <nbl>  - 


__£b 


5Z(5) 


)][, 


(B.10) 


where  § = cu  - kVz/ku^  and  Z(?)  is  the  usual  plasma  dispersion  func- 
tion8. The  important  result  here  is  that  ion  beam  contribution  to  the 
dispersion  equation  is  reduced  by  a factor  of  V2Qro/u^,  where 
V2/u^  « 1.  Terms  of  0(v*/u£)  or  smaller  may  be  ignored. 
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